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ABSTRACT 

The Subaru Deep Field (SDF) Supernova Survey discovered 10 Type la su- 
pernovae (SNe la) in the redshift range 1.5 < z < 2.0, as determined solely 
i- £h ■ from photometric redshifts of the host galaxies. However, photometric redshifts 

might be biased, and the SN sample could be contaminated by active galactic 
nuclei (AGNs). Unfortunately, measuring spectroscopic redshifts of galaxies in 
the "redshift desert" 1.5 < z < 2.0 is hard because any prominent emission lines 
get shifted out of the optical and into the near infrared. Here we report the first 
robust redshift measurement and classification of hSDF0705.25, an SDF SN la 
host galaxy. Using the X-shooter spectrograph on the Very Large Telescope, we 
measure a spectroscopic redshift of 1.5456 ±0.0003, consistent with its photomet- 
ric redshift of 1.552 ± 0.018. From the strong emission-line spectrum we are able 
to rule out AGN activity and show that the SN host galaxy is a low-metallicity, 
starbursting dwarf galaxy, similar to typical SN la hosts at lower redshifts. This 
represents a successful test case for the X-shooter spectrograph and shows that, 
though still hard, it is now possible to measure robust redshifts for SN host 
galaxies in the redshift desert. 
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INTRODUCTION 



The nat ure of the pr oge nitor stellar systems of ty pe la supernovae (SNe la) remains a 
mystery (see lHowellll201ll and lMaoz fc Mannuccill2012l for reviews). While bot h circumstan- 



tial and direct lines of e vidence point to a carbon -oxygen white dwarf (WD; iNugent et al. 



201 ll ; iBloom et al.ll2012l and see lLeibundgutll2000l for a review) as the progenitor, the other- 
wise stable WD must be ignited. The current consensus is that the carbon in the core of the 
WD is ignited due to the buildup of pressure, or temperature, resulting from mass accretion 
from a companion star in a binary system. The two leading scena rios for the nature of 



the progenito r binary system are the single degenerate scenario (SD; IWhelan &: Ibenlll973 



Nomotolll982l ). which contends that the WD accretes mass from a main-sequence, helium, o r 
giant star; and the double degenerate scenario (DP; llben fc Tutukovlll984l ; IWebbinklll984l ). 
in which the WD merges with a second CO WD through loss of angular momentum and 
energy to gravitational waves. 

One way to discriminate between the SD and DD scenarios is through the delay-time dis- 
tribution (DTD; the distribution of times that elapse between a short burst of star formation 
and the subseque nt SN la explosions). Whereas most SD DTDs fall to zero after a few Gyr 
(though see, e.g.. lDi Stefano. Voss. &: Claeysll201ll ; lHachisu. Kato. &: Nomotoll2012l ). the DD 
DTD is a power law with index ~ — 1, with power out to a Hubble Tim e. The SN la DTD 
has b een recovered dir ectly using different methods and SN samples (se e 



Maoz fc Mannucci 



20121 for a review, and lMaoz. Mannucci. fc Brandtll2012r iGraur fc Maozll2012l for subsequent 
studies). There are also indirect methods to recover the DTD. One such method relies on 
the measurement of volumetric SN la rates in field galaxies at cosmic time t, i?i a (t). These 
SN la rates are the convolution of the cosmic star formation history (SFH), S(t), with the 
DTD, 

Ru(t)= [ S(t-t')y(t')dt'. (l) 



Thus, the DTD can be recovered by comparing the field SN la rates to the cosmic SFH at 
different redshifts. 

Measurements of the volumetric SN la rate in field galaxies agree out to z ~ 1 [e.g., 
Perrett et al.ll2012t IGraur et al.ll201ll (Gil) and refer ences therein]. These measurements 
were first extended to z > 1 b y iDahlen et al.l (120041 ). with additional data analyzed by 



Dahlen. Strolger. fc Riessl (120081 ) . usi ng the Hubble Sp ace Telescope Advanced Camera for 
Surveys to survey the GOODS fields ( IRiess et al.l 120041 ) . Based on a sample of 3 SNe la at 
z > 1.4, they argued that th e SN l a rat e declined at z > 0.8. Fit t ing th is declining SN la 
rate evolution, IStrolger et al.l ( 120041 ) and lStrolger. Dahlen. fc Riessl (120 101 ) surmised that the 
DTD is confined to delay times of 3-4 Gyr. Conducting a similar SN survey in the Subaru 
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Deep Field (SDF) with the Subaru 8.2-m Telescope, Gil discovered 150 SNe, of which 28 
(10) were 1 < z < 1.5 {z > 1.5) SNe la. The Gil SN la rates indicate that the SN la rate 
evolution does not decline at high redshifts, but rather levels off, as predicted by the DD 
scenario. 

However, the classification of the SNe discovered in the SDF is purely photometric and 
depends on the redshift of the host galaxy. As detailed below, the redshifts of most of the 
SDF SN host galaxies, including those at z > 1.5, are photometric redshifts (photo-z's). 
These photo- z's are based on photometry in 11 bands, from the far-ultraviolet (UV) to 
the near-infrared (IR), trained on hundreds of galaxies in the field with spectroscopy. Yet, 
because of the inherent difficulty in obtaining spectroscopic redshifts for early-type galaxies 
due to the lack of strong emission lines, training the photo- z method used by Gil in that 
range is difficult. There could be systematic biases in the redshift estimates in that bin, 
biases that are not accounted for in the formal uncertainty. Because of the small number 
of SN candidates, even a few 'catastrophic' photo- z failures could strongly distort the final 
results. While Gil used several methods to weed out interloping active galactic nucelei 
(AGNs), as each SN in the SDF sample was only observed in one epoch, there could still 
be some AGN contamination as well. In order to determine whether the Gil z > 1.5 rate 
suffers from such systematic biases, the host galaxies of the SN candidates must be observed 
spectroscopically, in order to measure their redshifts robustly, and to ascertain whether any 
of them harbor an AGN. 

In this paper we present a VLT/X-shooter emission-line spectrum, derive the spec- 
troscopic redshift and classify AGN vs. star-formation activity of the SN host galaxy 
hSDF0705.25. We also place constraints on the metallicity and star-formation rate of the 
host galaxy. 

Throughout this paper we assume a flat ACDM cosmology with H = 70 km s" 1 and 
Vt m = 0.3. 



2. DATA 



2.1. Photometric Data 



Gil derived SN host-galaxy photo- z's using the Zurich Extragalactic Bayesian Redshift 
Analyzer (zebra§; Feldmann et al. 2006 ). Briefly, the host galaxies are first fit to the Bemtez 



2 http:/ /www.exp-astro.phys.ethz.ch/ZEBRA/ 
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( 120001 ) galaxy spectral energy distributions (based on galaxy spectra from Coleman, Wu, 
& Weedman 1996 for E/SO, Sbc, Scd, and Irr galaxies and Kinney et al. 1996 for two 
starbursting galaxies). The residuals of these fits are analyzed for systematic errors, which 
are then removed from the host galaxy photometry catalog. This catalog was composed 
of photometry in 11 bands, from the far-UV to the near-IR: far-U V (FUV) and near-UV 



(NUV) from the Galaxy Evolution Exvlorer ( G ALEX:\Lv et al-lbopgh: B, V, R, i', z 1 , NB816 



and NB921 from the Subaru 8.2-m Telescope ( iKashikawa et al.ll2004f); and J and K, where 



available, from the United Kingdom Infrared Telescope ( lHayashi et al.l l2009t Motohara et 
al., in preparation). Next, the basic galaxy templates are corrected using a training set 
of galaxies with known spectroscopic redshifts. Gil used a training set of 431 galaxies, of 
which 150 were in the range 1 < z < 2 (see their section 4.2 for a detailed description). The 
resulting corrected templates are interpolated to produce a range of intermediate galaxy 
spectral templates. These templates are then used in an iterative Bayesian procedure, where 
the redshift and template probability distributions (PDFs) produced in each step are used 
as priors for the next one. 

The resultant root-mean-square scatter between the spectroscopic and photometric red- 
shifts of the Gil training set was 0a*/(1 + z s ) = 0.075 (where Az = z s — z p ) in the range 
< z p < 2, after rejecting six 4cr outliers. For the 24 Gil host galaxies with spectroscopic 
redshifts, the scatter was smaller: <ta*/(1 + z s ) = 0.044, with no 4er outliers. However, due 
to the dearth of spectroscopic redshifts at z > 1.5 (only 24, or ~ 6% of the training-set 
galaxies), there is a possible bias towards systematically higher values in the Gil z > 1.5 
photo- z's. 

Using the above method, the i'= 24 mag host galaxy of SNSDF0705.25, denoted hSDF0705.25, 
was typed as a Sbc galaxy with a sharp z-PDF that peaked at z p = 1.552±0.018. In Figured] 
we show the photometry of hSDF0705.25, along with the best-fitting galaxy spectral-energy 
distribution and resultant z-PDF. 

Each SN candidate in the SDF was observed in one of four independent epochs, in the 
R, i', and z 1 bands. Consequently, variable AGNs can be mistaken for SNe. Gil identified 
interloping AGNs using a catalog of known variable AGNs in the SDF (provided by T. 
Morokuma) and by culling SN candidates that appeared in more than one of the four survey 
epochs (see their section 3.1 for a detailed description). However, lacking spectroscopy, the 
possibility that some of the z > 1.5 SDF SNe la are in fact AGNs cannot be ruled out. 
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hSDF0705.25: z p = 1.552 ± 0.018, \ 2 =1-4, M B = -20.71 




0.5 1 1.5 2 1 1.2 1.4 1.6 1.8 2 

A [/im] Redshift 



Fig. 1. — zebra fit and resultant redshift PDF of hSDF0705.25. The left panel shows the ac- 
tual photometry (filled circles), the best-fitting galaxy template (solid line), and its synthetic 
photometry (empty circles). The vertical error bars denote the photometric uncertainty, and 
the horizontal error bars show the width of the filter. The la upper limit on the photometry 
in the GALEX FUVhaxid is shown as the downturned arrow. The header gives the designa- 
tion of the SN host galaxy, most probable photo- z (z p ), the \ 2 P er degree of freedom of the 
fit, and the absolute i?-band magnitude the object would have at z p . The right panel shows 
the resultant z-PDF. The red line marks the spectroscopic redshift obtained in this work. 
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Fig. 2. — X-shooter ID and 2D spectrum of hSDF0705.25. The spectrum has been corrected 
for Galactic extinction. The solid (red) line shows the fit to the emission lines assuming a 
Gaussian line profile. The dashed line shows a scaled version of the Ha line assuming the 
same FWHM (in velocity units) and extinction in the host of E(B — V) = mag (upper) 
and EiB — V) = 0.5 mag (lower). Color version available online. 
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2.2. Spectroscopic data 



The s pectrum of hSDF0705 . 25 was obtained on 2012 April 23 with the X-shooter spec- 
trograph fjD'Odorico et al.ll2006t IVernet et al.ll201ll ) at the Very Large Telescope (VLT) at 
Cerro Paranal, Chile. We used an ABBA on-source nodding template with an exposure time 
of 1.3 hr (4x1200 sec) and a 0"9 slhfl. The spectrum was obtained unde r clear conditions 
For de tails on high-redshift emission-line spectroscopy with X-shooter see iFrederiksen et al. 
fl2012h . 



The X-shooter spectra were reduced using the official X-shooter pipelin^f] vl.3.7. The 
extraction of the object spectrum was conducted with our own IDL script, and flux calibra- 
tion was done using the flux standard star, LTT 3218^1 



3. ANALYSIS 



The flux-calibrated spectrum is corrected for Galactic extinction and slit loss. The 

Gala ctic extinction along the line of sight to hSDF07 05.25 is Ay = 0.042 mag^l ( Schlafly fc Finkbeiner 



20 111 ) and we assume the Galactic extinction law of iFitzpatrickl ( 119991 ). The slit loss correc- 
tion assumes a A -a2 variation in seeing. 

In the available atmospheric transmission windows, we detect a strong emission line 
at 1671.88 nm which we identify as Ha at z = 1.54563 ± 0.00027. This identification is 
supported by the detection of H/3, [N II] 6583, and [S II] 6716,6731 at the same redshift. 
We do not detect any continuum in the spectrum, only the strong emission lines in Figure 
|2j The detection of H/3 is severely affected by noise and therefore not fitted. We do not 
detect [O III] 4959,5007 in the spectrum and can therefore only derive upper limits on the 
flux. The detected emission lines in the corrected spectrum are fitted individually with a 
Gaussian line profile and the flux is calculated from the fit. The redshift is determined from 
the centroid of the fit to the individual lines (see Table [I]) . We apply a weighted average 
with weights equal to the inverse of the variance on the centroids, w = a~ 2 . The derived 
redshift is consistent with the Gil photo-z value of z = 1.552 ± 0.018. 



3 A l'.'O slit in the UVB arm. 



Sec, http : //www . eso . org/ sci/ software/pipelines/ 



5 See http : //www . eso . org/sci/observing/tools/ standards/ spectra/ltt3218 .html 

6 Quoted from the NASA/IPAC Extragalactic Database (NED) 
http : //ned . ipac . caltech . edu/ 



website: 
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We plot in Figure[2]a scaled-down version of the Ha line assuming case-B recombination, 
the same F WHM (in veloc ity units) and intrinsic reddening in the host of E(B — V) of zero 
or 0.5 mag ( jCafzettill200ll ) to illustrate the variation allowed by the spectrum. 



The emission line ratios N2 = log([N II] 6583/Ha) and 03 = log([0 III] 5007/H/3) are 
among the main diagnostics for discrimina t ing; b etween star formation and AGN activity by 
way of the Baldwin. Phillips. &: Terlevichl ( 119811 . BPT) diagram. The upper limit on 03 is 
calculated using the 3a upper limit on [0 III] and the predicted H/3 flux from the Ha flux. 
We calculate the H/3 flux assuming E(B — V) = or 0.5 mag to make sure the classification 
does not depend on the assumed amount of extinction in the host galaxy. In both cases the 
host galaxy falls within the star- forming region of the BPT diagram (see iFrederiksen et al. 
20121 . their Figure 5). The lack of a broad line component in the emission lines (see Tabled]) 
also supports the conclusion that the line-flux is dominated by star formation and not AGN 
activity. 

To constrain the metallicity of the host galaxy we make use of the emission-line ra- 
tios N2 and iV203 = lo g([N II] 6583/ [O III] 5007). Using the emission-line calibrations of 
Kewlev fc Dopital d2002h . we apply a three-step iteration to set an upper limit on the metal- 
licity. The first step involves using the N2 ratio to constrain the metallicity to either the 
upper or lower metallicity regime. The N2 ratio has the advantage of being virtually insen- 
sitive to extinction due to the close proximity of the two emission lines. The N2 ratio places 
the host galaxy in the low-met allicity regime at 12 + \og{0/H) < 8.55. Next, we use the 
iV203 ratio to refine the initial estimate. An upper limit on the [O III] line flux translates 
into a lower limit on iV203. The lower limit on iV203 places an upper limit on the ionization 
parameter q < 2 x 10 7 cm s _1 , but due to the uncertainty in the intrinsic extinction of the 
host we take the precaution of placing the upper limit at q < 4 x 10 7 cm s" 1 , which would 
be consistent with EiB — V) ~ 0.25. Finally, we refine the metallicity estimate using the 
constraint on q. This gives us our final estimate of 12 + \og(Q/H) < 8.0. For compariso n, 



the metallicity derived from the [N II]/Ha line ratio is 8.40 ± 0.18 (IPettini fc Page! 



200J) 



We convert the Ha flux to a res t-frame Ha lumin osity and furth er into a s t ar-for mation 
rate (SFR) using the calibration of iKennicuttl (119981 ). rescaled to a IChabrierl (120031 ) initial 
mass function (IMF). The observed SFR, 4.0 ± 0.7 M yr _1 , represents a lower limit to the 
intrinsic SFR of the host galaxy, due to the unconstrained extinction in the host galaxy. 

The Photo-z fitter zebra is not suited for deriving stellar parameters like the mass of the 
galaxy. We therefore fit the Subaru+UKIRT (B, V, R, i', 2!, J and K) photometric points 
(corre cted for foreground extinction in the same way as the spectrum) using the FAST SED 
fitter (IKriek et al.ll2009i ). We derive the intrinsic extinction, stellar mass and stellar age 
of the host galaxy (see Table [2]). The extinction in the host is not very well constrained 
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(0 < Ay < 1.1) so we assume Ay = in our further analysis. We calculate the specific SFR, 
sSFR, using the Ha SFR and the stellar mass from FAST. 



DISCUSSION AND CONCLUSIONS 



Our derived spectroscopic redshift, z = 1.54563 ± 0.00027, is in full agreement with the 
Gil photometric redshift of z = 1.552 ± 0.018. We also exclude AGN activity as the source 
of the emission-line flux. From the flux of the Ha line we derive an observed SFR and from 
the emission line ratios we constrain the host-galaxy metallicity. 

We find a highly star-forming (i.e., high sSFR), low- metallicity SN host galaxy (see 
Frederiksen et al.l 12012k for the treatment of another high SF low metallicity SN host at 
similar redshift). The derived SFR is affected by extinction in the host galaxy, but due to 
the low signal-to-noise ratio of the H/3 line we are not able to place any strong constraints 
on it. Likewise the extinction derived from SED fitting does not provide a strong constraint. 

The stellar mass and SFR place hSDF0705 .25 midway betwee n th e z = 1 and z = 2 



Daddi et al. 


(2007 


) and 


Elbaz et al. 


(2007) 



As hSDF0705.25 falls within the main sequence at its redshift, it can be classified as an 
average star- forming galaxy at its redshift. Alternatively the specific SFR (sSFR) can be 
used to define w hether a galaxy is a passive, a star forming or a star burst galaxy. Using 



the definition of ISullivan et al.l ( 120061 . 0.2 < z < 0.75, see their Figure 6), the high sSFR of 



hSDF0705.25 makes it a star-burst galaxy. 

Having two different measures of metallicity that do not agree is not unusual, as the 
absol ute value of the me t allicity depends on the calibration (also called the metallicity scale) 
used. iKewley fc Ellison! (120081 ) made the comparison between different metallicity calibra- 
tions and determined transformations to translate from one metallicity scale to another. 
From the transformation we expect PP04 to be systematically lower then KD02. The metal- 
licity derived for this host galaxy is inconsistent with this tendency as the upper limit on the 
KD02 metallicity is ~ 2a lower then the PP04 metallicity. This is can be due to extinction 
effects or to some differences between theoretical and empirical calibrations that are not 
captured in KE04. The extinction effect is most important for the N203 line ratio (due to 
the wide separation between the lines). We have tried to account for extinction by assuming 
a conservative upper limit on N203 (assuming E(B — V) — 0.25). We therefore believe that 
the difference is due to differences between theoretical and empirical calibrations. 

SNe la from low mass galaxies w i th high sSF R (and therefore low metallicity, see 
Lara-Lopez et al.l |2010| ; iMannucci et al.l l2010l 1201 ll ) constitute the the bulk of the z < 1 
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Table 1. Spectroscopic summary 



Line 


Wavelength 


Obs. FWHM a 


Int. FWHM a 


Flux 




nm 


km s _1 


km s" 1 


cr g g-1 cm -2 


Ha 


1671.88 ±0.09 


117 ±18 


103 ±18 


58.8 ± 10.4 


[N II] 6583 


1677.45 ± 0.07 


47 ±14 




7.7 ± 3.1 


[S II] 6716 


1711.04 ±0.11 


50 ±21 




10.4 ± 5.0 


[S II] 6731 


1714.64 ±0.09 


73 ±16 


47 ± 16 


9.9 ± 2.8 


[0 III] 5007 b 








< 23.2 


[0 III] 4959 b 








< 9.1 



a Thc instrumental resolution is 57 km s 
b This is the 3a upper limit on the flux. 



Table 2. Derived properties of SN host galaxy hSDF0705.25 



Property Value 

Redshift z = 1.5456 ± 0.0003 

Star-formation rate ab SFR = 4.0 ± 0.7 M Q yr" 1 

Metallicity (PP04) 12 + log(0/H) = 8.4 ± 0.2 

Metallicity (KD02) 12 + log(0/H) < 8.0 

Ionization parameter (KD02) q < 4 x 10 7 cm s -1 

Stellar mass c log(M„[M Q ]) = 9.46±g;g| 

Specific SFR log^SF^Kyr- 1 ]) = -8.86^ 

c A___ _ n /I +0.7 



Host extinction c v4yhost = 0.4_ ; 4 mag 

h0.3 
-0.6 



Stellar age c log(t*[yr]) = 7.9 +0 ' 3 



a Assuming a Chabrier IMF. 

b SFR is a lower limit as extinction in the host galaxy would make the intrinsic SFR higher. 

c Using an exponentially declining star-formation history, redshift fixed to the spectroscopic redshift, and 
metallicity fixed to Z = 0.008 (12 + log(0/H) - 8.4). 
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SN la population (jSullivan et al.ll2006l see). hSDF0705.25 is a similar low-metallicity, high- 
sSFR, host galaxy at z — 1.55 and therefore indicates the presence of such "prompt" SNe la 
also at these higher redshifts. 

The SN la rates at high redshift are dominated by small-number statistics. The three 
z > 1.4 HST/GOODS SNe la belong to host galaxies with measured spectroscopic redshifts 
and no AGN activity. On the other hand, the larger Subaru/SDF sample includes ten 
SNe la at z > 1.5, but their classification as SNe la relies on photometric redshifts which, 
at high redshifts, might be systematically offset. Two ongoing HST Mult i- Cycle Treasury 
programs, the Cluster Lensing and Supernova Survey with Hubble flPostman et al.l 1 2 1 21) 



and the Cosmic Assem bly Near-IR Deep Extragalactic L egacy Survey (iGrogin et a 



Koekemoer et al.l 1201 ll ) will find new SNe out to z ~ 2.5 flRodney et al.l 



2012 



Jones 



2011 



20121) 



but their samples will still be small and will suffer from the same classification challenges 
faced by the GOODS and SDF surveys. It is thus imperative to test the robustness of the 
SDF SN la rate measurements b y spectroscopically measuring the redshifts of the SN host 
galaxies (IFrederiksen et al.ll2012l ). and ascertaining whether the SN sample was contaminated 
by unclassified AGNs. The confirmation of the photometric redshift of hSDF0705.25, and 
its classification as a star-forming, non-AGN-hosting galaxy is the first step along this road. 

Galaxy redshifts are usually measured either from bright emission lines or from the 
4,000 A break. Beyond redshift z ~ 1.5, the strong emission lines, like H/3 and [O III] 4959+ 
5007, that we rely on to measure a galaxy's redshift are shifted into the near-IR and out of the 
wavelength range covered by most existing spectrographs, leaving only a single [O III] 3727 
line in the optical window, with no other strong emission lines to confirm the redshift. The 
signal-to-noise of the [O III] 3727 line depends on the metallicity of the object, the amount 
of telluric absorption, and the atmospheric airglow lines present in the red part of the optical 
windows. Determining the redshift from the continuum requires prohibitively long exposure 
times. This is why z ~ 1.5 marks the start of the "redshift desert," where to date it has 
been exceedingly hard to measure galaxy redshifts. We have shown that X-shooter, with its 
uniquely long wavelength coverage from the near-UV to the near-IR, makes i t possible to 
determine robust redshift measurements well into the "redshift desert" (see also Christensen 



20 111 ; iKriihler et al.ll2012l ). Yet this is still far from an easy task. hSDF0705.25 has been a 
test case and the first of the ten z > 1.5 SDF SN la host galaxies to receive a spectroscopic 
redshift. Of the remaining nine galaxies, three more can be targeted by X-shooter, two are 
too faint, and four are passive galaxies with no predicted emission lines. For these galaxies, 
the 4,000 A break can only be targeted with the sensitivity of the Hubble Space Telescope. 



We thank Martin Sparre for providing his X-shooter meta-pipeline, which has simplified 
the reduction of the X-shooter spectra significantly. The Dark Cosmology Centre is funded by 
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